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ABSTRACT 

A two-dimensional computer code to solve the Bur- 
nett equations has been developed which computes the 
flow interaction between an exhausted plume and hy- 
personic external flow near the afterbody of a flight ve- 
hicle. This Burnett-2D code extends the capability of 
Navier-Stokes solver (RPLUS2D code) to include high- 
order Burnett source terms and slip-wall conditions for 
velocity and temperature. Higher-order Burnett viscous 
stress and heat flux terms are discretized using central- 
differencing and treated as source terms. Blocking logic 
is adopted in order to overcome the difficulty of grid gen- 
eration. The computation of exhaust plume flow field is 
divided into two steps. In the first step, the thruster noz- 
zle exit conditions are computed which generates inflow 
conditions in the base area near the afterbody. Results 
demonstrated that at high altitudes, the computations 
of nozzle exit conditions must include the effects of base 
flow since significant expansion exists in the base region. 
In the second step, Burnett equations were solved for 
exhaust plume flow field near the afterbody. The free 
stream conditions are set at an altitude equal to 80 km 
and the Mach number is equal to 5.0. The preliminary 
results show that the plume expansion, as altitudes in- 
creases, will eventually cause upstream flow separation. 

INTRODUCTION 

For a single-stage-to-orbit hypersonic flight vehicle 
cruising at high altitudes, an exhaust plume is created 
in the base area by a propulsive device. This plume in- 
teracts with the free stream near the afterbody and may 
create the undesirable plume-induced separation. This 
separation affects the dynamic characteristics of the flight 
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vehicle. Due to the difficulty of collecting reliable wind- 
tunnel data, numerically simulated data has become im- 
perative for this purpose. 

As the flight altitude increases from sea level to 
about 70 km, the plume/free-stream interaction flow field 
is in the continuum regime. The conventional Navier- 
Stokes equations are valid governing equations. As al- 
titude increases from 70 to 100 km, the ambient air 
becomes more rarefied and the flow is usually in the 
transitional regime. The transitional flow regime is de- 
fined when the Knudsen number Kn is from 0.01 to 
1.0. The prevailing reason for the failure of applying the 
Navier-Stokes equations in the transitional flow regime 
is primarily due to the breakdown of the translational 
equilibrium^ 1,2 !. it is common to apply the Direct Sim- 
ulation Monte Carlo (DSMC) method for rarefied gas- 
dynamic problems. Moss^ 3 ! conducted a series of DSMC 
analysis and concluded that the Navier-Stokes equations 
cease to be valid above 92 kilometers. The drawback 
of the DSMC calculation for flow problems with small 
Knudsen number is that the computational time becomes 
prohibitively expensive. 

To circumvent these discrepancies and to preserve 
the continuum modeling, the third-order Burnett equa- 
tions are hypothesized as the governing equations in the 
transitional regime. Zhong^ 4,5 ! conducted a series of hy- 
personic flow analysis for flow over a blunt body us- 
ing the Augmented Burnett equations. Liaw, Deng and 
Chout 6,7,8 ! a lso conducted a series of analysis for shock 
wave structure prediction and low-density nozzle flow cal- 
culations to solve the Burnett and Boltzmann equations. 
It was concluded that the Burnett solutions are more 
accurate than the Navier-Stokes solutions in the transi- 
tional regime. 

The objective of the current research is to numeri- 
cally investigate the exhausted plume flow field near the 
afterbody of a hypersonic flight vehicle in high altitudes 
by solving the Burnett equations. 

MATHEMATICAL FORMULATIONS 

The Burnett equations are derived to include the 
second-order term of the distribution function in the 
Chapman-Enskog expansion series, 
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/ = fo + A'n/i + Kn 2 f 2 + • * *, (1) 

so that the Knudsen number effect is emphatically em- 
bedded into the moment equations. As far as the math- 
ematical manipulation is concerned, the Navier-Stokes 
equations are the subset of the Burnett equations be- 
cause the Navier-Stokes equations are obtained with the 
first order Chapman-Enskog expansion. In consequence, 
the Burnett equations are a system of third-order non- 
linear partial differential equations which are far more 
complicated than the Navier-Stokes equations. 

Written in a general body-fitted coordinate system, 
(£, r\) , the two-dimensional Burnett equations can be ex- 
pressed as 
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coordinates (x,y), p is pressure, e t is the specific total 
energy, r xx , r xy , r yy are the Navier-Stokes shear stress 
components, 7xi\ r xy \ r yy ^ are the Burnett stress com- 
ponents, q Xy q y are the first-order heat flux components, 
and qi 2 \ <j y 2) are the Burnett heat flux components. 

The Burnett stress components (r xx , r xy , r yx , r yy ), 

/ o \ ( o \ 

and heat flux components (q x \q y ), are derived from 
the previous work^l. In two-dimensional Cartesian coor- 
dinates (x, y), they can be expressed as 
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In equations (3), (4), (5), (6) and (7), H is the cell volume, 
p is density, u and v are velocity components in Cartesian 
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where, the subindex x and y represent the partial deriva- 
tives with respect to coordinates x and y, e.g., 

8T dT _ d 2 T 
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The coefficients 7 * are obtained based on the 

Chapman-Enskog expansion and are depend on the 
molecular force model^. In this research, the 
Maxwellian gas model is assumed. 


NUMERICAL PROCEDURES 


Numerical Scheme for Burnett Equations. The 
solution methodology for two-dimensional Burnett equa- 
tions is extended from the RPLUS2D^ code. The 
second-order upwind (UP) scheme, namely, the Steger- 
Warming vector splitting scheme, and the Lower-Upper 
(LU) approximate factorization methods are used to dis- 
cretize the Burnett equations. The discretized Burnett 
equations can be written as 

(I + Ar(A-i+ J+ A+i- + A-J3+. + A+S-.)6Q 

= —A rResi j 

(13) 


where, the right-hand-side residue Resij is 

Re,, i =^(E-E.) + ±(F-h) (14) 


In equations (13) and (14), A+ and A are forward and 
backward one-sided differencing operators, and 
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in which A^,Aj}V * * etc. are the diagonal eigen-value 
matrices, and etc. are the right eigen-matrices ob- 
tained from the Steger- Warming splitting. Substituting 
A+ and A“ into equation (14), it becomes 

^D + L + U^j 6 Q = -AtResij (16) 
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After applying approximate factorization procedures, it 
becomes 
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where D is a diagonal matrix, L and U are upper and 
lower matrices, respectively. At each time step, the op- 
eration of LU-SW/UP scheme can be written as: 

(D + L) 6 Qlj = — At Re Si j 

(D + U) 6 Qij = DSQij (18) 
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For residual ( Resi j ) calculations, the Osher s up- 
wind differencing scheme is used for inviscid terms, and 
central differencing is used for viscous and Burnett terms. 
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The Burnett stress and heat flux terms are treated 
as source terms in the solution procedure. At each time 
step, the source terms are calculated explicitly based 
on flow parameters of the previous time step. A time- 
marching technique is used to obtain the steady state 
solutions. A steady state solution is obtained when all 
residues (density, velocity and enthalpy) are less than 
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Blocking Algorithm and Boundary Conditions. 
To reduce the difficulty of grid generation, the blocking 
logic was used to compute the flow near the afterbody. 
The computational domain is divided into two blocks of 
grid system as shown in Figure 1. The first block is con- 
structed in the lower part of the afterbody. The second 
block is constructed in the plume region near the after- 
body. For block-one, free stream, wall, interface and su- 
personic outflow boundary conditions must be specified 
for equilibrium air flow. For block-two, thruster noz- 
zle exit conditions, wall, interface and outflow conditions 
need to be given. The thruster nozzle exit conditions 
were computed independently using the RPLUS2D code. 
The nozzle calculation is assumed a H 2 /O 2 equilibrium 
flow. The base flow area is included in the computational 
domain and the flow field in the base area is computed 
simutaneously with the nozzle internal flow field. The 
computed nozzle exit conditions are then averaged and 
used as the inflow boundary for block-two. The inter- 
face conditions have to be coupled between block-one and 
block-two. A first-order upwind approximation scheme 
is applied at the interface. Since the conventional con- 
tinuum no-slip wall conditions break down as the flow 
becomes rarefied, the slip- wall conditions are applied at 
the wall. Based on the first-order Maxwell/Smoluchowski 
slip boundary conditions, the wall temperature T s and 
slip velocity u, (Schaaft 10 !) are derived as 
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where 5 is the surface tangential direction, n is the inward 
surface normal direction, T w is the wall temperature, and 
<r n ,cr f are the surface reflection and accommodation fac- 
tors, respectively. A complete accommodation, <r n = 1, 
<r t = 1, is assumed in this research. The slip wall con- 
ditions are implemented at each surface point implicitly. 
The calculated slip velocity and surface temperatures are 
then utilized to compute boundary fluxes. The velocity 
and temperature in the boundary cells are calculated as 

boundary cell = 2 V$lip V| nternal cell (22) 

T boundary cell = 2 T s — T{ n ternal cell 


for ambient pressures equal 14.7, 1.0 and Chips*, respec- 
tively. At sea level, the typical plume expansion is shown 
in Figure 2(a). As altitude increases, the ambient pres- 
sure decreases. At ambient pressure equals l.Opst, ex- 
haust plume expands near the nozzle lip. Shock waves 
are formed near the edge of the plume. Further increas- 
ing altitude will generate higher expansion plume in the 
base area of the nozzle. Physically, the high expansion 
plume and the boundary layer near the nozzle lip must 
be calculate simultaneously during numerical iterations. 
However, it was neglected for simplicity. These results 
demonstrated that in order to obtain accurate nozzle exit 
conditions, base area and external free stream must be 
computed simultaneously with the nozzle internal flow 
field. 


Since the solution methodology of BURNETT-2D 
is adopted from the RPLUS-2D code, detailed descrip- 
tions of the numerical algorithm can be referenced in 
Shuent 11,12 ’ 13 !, and Tsai and Hsieh^. 

RESULTS AND DISCUSSION 

To avoid computational complexity, the exhaust 
plume flow field near the afterbody of a hypersonic flight 
vehicle at a high altitude is computed through two steps. 
First, the thruster nozzle exit conditions are calculated 
through the Navier-Stokes solver - RPLUS2D. Secondly, 
the Burnett equations are solved for exhaust plume flow 
field near the afterbody. The computed nozzle exit condi- 
tions from the first step are averaged and used as input 
conditions for free stream/plume interaction. This av- 
eraging process has simplified the global computational 
domain, but the boundary layer effects are neglected. 

Nozzle Exit Conditions Computation. The thruster 
nozzle has a throat diameter of 6 inches with an area 
ratio of 11.4. The chamber pressure is 1600psi and the 
chamber temperature is 4650/?. An 18 step hydrogen/air 
chemistry model is applied to compute the nozzle inter- 
nal flow field. The gas mixture is assumed to have 9 
species. As altitude increases above 80A;m, the air be- 
comes highly rarefied. Rapid expansion may occur in 
the base region. This highly expanded plume will inter- 
act with free stream and the flow in the forward region, 
which in turn, affect the nozzle exit conditions. In order 
to provide realistic boundary conditions at nozzle exit for 
the afterbody flow field calculation, flow fields in the base 
area and region outside the nozzle are computed simul- 
taneously. The nozzle wall lip effects must be also taken 
into account in order to calculate realistic nozzle exit 
conditions. For preliminary study, the nozzle lip effect 
is neglected and the ambient free stream Mach number 
is set equal to 3.0 for all cases. It was assumed that no 
chemical reaction occurs in the base area and outside re- 
gion of the nozzle, and the gas mixture contains only air. 
RPLUS2D code was used to compute the entire flow field. 
Three ambient pressures were selected for the calculation. 
Figures 2(a), 2(b) and 2(c) show the Mach number con- 
tours of the nozzle plume and free stream interactions 


Plume Flow Field Prediction near the Aftbody 
Solving the Burnett Equations. A hypothetical two- 
dimensional body is shown in Figure 3(a). A high ex- 
pansion ratio thruster nozzle is placed inside the lower 
part of the body and used as a propulsive device. The 
computed nozzle exit conditions from step one are aver- 
aged and used as input conditions at the inflow boundary. 
Since the ambient conditions for the nozzle calculations 
are unknown, the results from step one can only be used 
as a reference. Further study is needed in order to obtain 
better nozzle exit conditions. Three cases were computed 
for exhaust plume flow field near the afterbody. The 
nozzle exit static pressures are specified as 5, 10 and 20 
times of the free stream value. Based on results from 
step one, it is estimated that the nozzle exit Mach num- 
ber at an altitude of 80ibm is about 3.0 for the selected 
nozzle configuration. Figures 3(a), 3(b), 3(c) show the 
two-dimensional Burnett solutions of exhaust plume flow 
field near the afterbody of a two-dimensional flight vehi- 
cle. The free stream Mach number is 5.0 and the flight 
altitude is 80Arm for all three cases. A strong oblique 
shock is formed near the leading edge of the flight vehi- 
cle. Exhaust plume from nozzle exit interacts with the 
wall and free stream to form a shockwave/viscous inter- 
action flow pattern near the afterbody. With increasing 
nozzle exit pressure, the exhaust plume expands toward 
the upstream. It is expected that further increasing of 
the nozzle exit pressure will introduce a plume-induced 
separation near the afterbody. 

CONCLUSION 

A two-dimensional code to solve the Burnett equa- 
tions for free stream/plume interaction near the after- 
body of a hypersonic flight vehicle was developed. The 
commutation of exhaust plume flow field near the after- 
body contains two steps. The thruster nozzle exit condi- 
tions were determined a priori through a Navier-Stokes 
nozzle solutions which provide upstream boundary con- 
ditions for the global flow field calculation of the free 
stream/plume interaction. Results demonstrated that at 
high altitudes, the base region and outside region must 
be included for calculations of nozzle exit conditions si- 



multaneously. Exhaust plume flow field at an altitude 
of 80fcm near the afterbody was successfully computed 
using the developed Burnett code. Results showed that 
with increasing nozzle exit pressure, the plume induced 
separation may occur near the afterbody which will sig- 
nificantly affect the aerodynamic characteristics of a hy- 
personic flight vehicle. 
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Figure 1. Schematic showing of grid setup for computation of exhaust plume flow field near the afterbody. 
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Figure 2(a). Ambient Pressure 14.7pst\ 
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Figure 2(b). Ambient Pressure l.Opsi. 
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Figure 2(c). Ambient Pressure O.lpsi. 

Figure 2. Mach number contours for the nozzle plume and free stream interaction. 
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Figure 3(b). Pressure ratio P e /Poo = 10. 
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Figure 3(c). Pressure ratio Pe/Poo = 20. 

Figure 3. Two-dimensional Burnett solutions of exhaust plume flow field near the aft body. 






